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Abstract. We present a detailed study of a statistically 
complete sample of highly inclined disk galaxies in the 
near-infrared K' band. Since the if'-band light is rela- 
tively insensitive to contamination by galactic dust, we 
have been able to follow the vertical light distributions all 
the way down to the galaxy planes. 

The mean levels for the sharpness of the X'-band lumi- 
nosity peaks indicate that the vertical luminosity distri- 
butions are more peaked than expected for the interme- 
diate sech{z) distribution, but rounder than exponential. 
After fitting a generalized family of fitting functions char- 
acterised by an exponent 2/n (n = oo for exponential, 
n = 2 for sech and n = 1 for sech^; van der Kruit 1988) 
we find that the mean value for 2/n in the K' band equals 
<2/n>K'— 0.538, CTif' — 0.198. Since projection of not 
completely edge-on galaxies onto the plane of the sky 
causes vertical luminosity profiles to become rounder, we 
have performed simulations that show that it is possible 
that all our galaxies can have intrinsically exponential ver- 
tical surface brightness distributions. 
We find that the profile shape is independent of galaxy 
type, and varies little with position along the major axis. 
The fact that we observe this in all our sample galaxies 
indicates that the formation process of the galaxy disks 
perpendicular to the galaxy planes is a process intrinsic 
to the disks themselves. 

Key words: galaxies: fundamental parameters — galax- 
ies: photometry — galaxies: spiral — galaxies: statistics 
— galaxies: structure 
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1.1. Motivation 

It is well-known from optical observations that many disk 
galaxies show a universal, exponential behaviour of the 
vertical luminosity away from the galaxy planes. This ob- 
servation, and the assumption that the mass-to-light ratio 
in the disk is constant, leads us to believe that galaxies 
have a universal, exponential vertical mass density profile 
(e.g., Tsikoudi 1979; van der Kruit & Searle 1981a). 

This nearly universal vertical structure is likely the re- 
sult of secular, internal evolution (van der Kruit & Searle 
1981a; Carlberg 1987). As Carlberg (1987) states, the con- 
tinuous variation of stellar kinematics from the youngest 
to the oldest disk stars strongly suggests that an ongoing 
dynamical evolution is indeed present. 

We know from observations of our Galaxy that galaxy 
disks probably consist of multiple components of increas- 
ing age, velocity dispersion and scale height; the den- 
sity is the sum of all these components. An exponen- 
tially decreasing density distribution with distance from 
the galaxy plane therefore puts interesting constraints on 
the star formation rate (SFR) and the dynamical evolu- 
tion of the disk. Burkert & Yoshii (1996) have shown, that 
these exponential vertical density profiles are a natural 
result of disk evolution if gaseous protodisks settle into 
isothermal equilibrium prior to star formation. If the star 
formation and cooling rates are comparable, stellar expo- 
nential 2-profiles arise due to the gravitational contraction 
of the gas towards the galaxy plane. 

1.2. Observational status 

Several models have been proposed and tested to account 
for the vertical distributions observed close to the galaxy 
planes. 

Van der Kruit & Searle (1981a,b, 1982a,b) studied the 
surface brightness distributions of edge-on disk galaxies 
in optical passbands, which were significantly affected by 
dust contamination. From these studies, they proposed 
the self-gravitating isothermal sheet (Spitzer 1942) as a 
model for the description of the vertical light distribution: 



L{z) = Lq sech^(z/zo) 



(1) 
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where Lq is the surface luminosity in the plane of the 
galaxy, zq is the vertical scale parameter and z is the 
distance from the galaxy plane, respectively. At large 
z heights, the vertical scale heights are equal for both 
the isothermal model and the exponential approximation, 
where zq equals twice the exponential scale height. At in- 
termediate z, most galaxies are dominated by the thin disk 
luminosity, which is thought to be locally isothermal (van 
der Kruit & Searle 1981a). 

The relations among the values for the density in the 
plane, p{0), the vertical scale parameter, zq, and the verti- 
cal velocity dispersion {v^)^^'^ for the solar neighbourhood 
show that the disk of our Galaxy is self-gravitating (van 
der Kruit & Searle 1981a). Based on the quite similar val- 
ues for Zq found by Van der Kruit & Searle (1981a) they 
argued that it may be reasonable to assume that the disks 
of other galaxies are also self-gravitating. 

Near-infrared observations of edge-on disk galaxies 
have shown an excess of light over the isothermal model at 
small distances from the galaxy planes, where the optical 
photometry is strongly affected by dust absorption (e.g., 
Wainscoat et al. 1989; Aoki et al. 1991; van Dokkum et al. 
1994). Wainscoat et al. (1989) show that the z-dependence 
of the light in the large southern edge-on IC 2531 demon- 
strates a more strongly peaked profile than expected from 
the isothermal sheet model, which appears to be better 
fitted by an exponential: 

L{z) = Lq exp{-z/h;,), (2) 

(where hz is the exponential vertical scale height), al- 
though their limited resolution does not unambiguously 
differentiate between the models. For the vertical iiT-band 
light distribution in NGC 891, Aoki et al. (1991) find that 
the exponential model fits the data remarkably well up to 
those 2-distances where the seeing convolution becomes 
significant. 

1.3. An intermediate solution 

Although the exponential model is mathematically attrac- 
tive because of its simplicity, there is no firm physical basis 
for such a model. An exponential vertical mass density dis- 
tribution can be constructed by adding up multiple stellar 
disk components. This can only be done if the contribu- 
tions from stars with larger velocity dispersions are in- 
creasingly dominating with increasing distance from the 
galaxy plane. However, a mechanism to account for such 
a process is as yet unknown (Burkert & Yoshii 1996). 

As van der Kruit (1988) argues, a pure exponential 
distribution also has some undesired properties, the most 
important one being a sharp minimum of the velocity dis- 
persion in the plane. Fuchs & Wielen's (1987) results show 
moderate gradients, much smaller than required for the 
exponential distribution (Bahcall 1984a,b). Therefore, van 
der Kruit (1988) proposed that an intermediate distribu- 



tion, such as the "sech(z)" distribution, could be a more 
appropriate one to use: 

L{z) = Lo sech{z/hz), (3) 

to account for the deviations from an isothermal sheet in 
the galaxy planes. 

Wainscoat et al.'s (1989) near-infrared photometry of 
IC 2531 would agree with this model, as would star counts 
for our Galaxy (e.g., Gilmore & Reid 1983). 

1.4. The Galaxy 

Studies of the Galaxy provide valuable information on the 
vertical structure of galaxy disks. These studies, based on 
star counts, have the advantage over the studies of exter- 
nal galaxies that they are less affected by dust absorption 
(in the solar neighbourhood) and effects of the presence of 
a young stellar population. Moreover, studies of the verti- 
cal luminosity structure in our Galaxy benefit greatly from 
the higher spatial resolution compared to that in external 
galaxies. Gilmore & Reid (1983) and Pritchet (1983) con- 
clude that the stellar z-distribution in our Galaxy is better 
approximated by an exponential rather than an isother- 
mal profile. 

On the other hand. Hill et al. (1979) derived density 
laws for A and F dwarfs towards the North Galactic Pole, 
which cannot be fit well by an exponential distribution, 
although they may be approximated as such in short dis- 
tance bins. Although they find that the F stars are roughly 
consistent with a single exponential, the A stars can only 
be approximated by a single exponential closer to the 
Galactic Plane and with a significantly smaller scale height 
than the F dwarfs. 

Based on observations in the near-infrared, Kent et al. 
(1991) concluded that the vertical light distribution (and 
hence probably also the mass distribution) follows an ex- 
ponential law more closely than an isothermal sheet ap- 
proximation. However, they did not compare the observed 
light distribution to other, intermediate models. 

1.5. Near-Infrared Observations 

The study of edge-on galaxies in the near-infrared is valu- 
able to reveal the true stellar distributions, as the near- 
infrared wavelengths permit to study these even at small 
z. 

In this paper we study the vertical luminosity profiles 
of a statistically complete sample of edge-on disk galaxies 
in the K' band. This sample is the largest of its kind 
available at the moment. For that reason we are able to 
study the structure of edge-on galaxies in a statistically 
consistent way, even down to very small distances from 
the galaxy planes. 

The JC'-band light (which is comparable to the stan- 
dard Johnson K band) is likely to be dominated by the old 
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disk, and the young disk contribution is relatively unim- 
portant. Moreover, the mass-to-light ratio is almost inde- 
pendent of metallicity and age in this band. Rix & Rieke 
(1993) find, from monitoring the gravity-sensitive CO(2.3 
/xm) index in the large disk galaxy M51, that young red 
supergiants do not distort the K-hand image significantly. 
They find that at most small portions of the spiral arms 
at have large contributions from young stars. 

The JT-band wavelength is too short for a substantial 
amount of direct emission from the dust. The high dust 
temperatures required to emit in K (800 - 1000 °K) are 
associated only with young stellar objects and compact 
Hii regions (see, e.g., Wainscoat et al. 1989). 

The light in Kis dominated by giants, which constitute 
only a small fraction of the stellar mass. However, old 
population giant stars have the same spatial distribution 
as the main-sequence stars (Rix & Rieke 1993). 

Therefore, since neither dust nor young, luminous red 
stars strongly affect the iiT-band image, i^'-band imaging 
with infrared arrays is a reliable and efficient method to 
map surface mass variations through surface brightness 
variations (Rix & Rieke 1993) 

In Sect. 1^ we present the sample properties and de- 
scribe the data reduction method used. The results of our 
detailed analysis of the vertical profiles are presented in 
Sect. ^ in which we also compare the results obtained in 
the optical / band to those from the K' band observations. 
In Sect. ^ we discuss the main observational results in the 
context of star formation and global galaxy structure pa- 
rameters. Finally, we summarize and conclude the paper 
in Sect. ||. 

2. Approach 

2.1. Sample selection, observations and data reduction 

To study the structural parameters of edge-on spiral galax- 
ies we selected a statistically complete sample taken from 
the Surface Photometry Catalogue of the ESO-Uppsala 
Galaxies (ESO-LV; Lauberts & Valentijn 1989) with the 
following properties: 

— their inclinations are greater than or equal to 87°; 

— the angular blue diameters (-D25) are larger than 2. '2; 

— the galaxy types range from SO to Sd, and 

— they should be non-interacting. 

The inclinations were determined following Guthrie 
(1992), assuming a true axial ratio logi?o = 0.95, cor- 
responding to an intrinsic flattening qo ~ (b/a)o of 0.11. 
From this intrinsic flattening the inclinations i were de- 
rived by using Hubble's (1926) formula 

COs'i = {q'-ql)/{l-q',), (4) 

where q = b/a is the observed axis ratio. 

Of the total sample of 93 southern edge-ons, an arbi- 
trary subsample of 24 galaxies was observed in the near- 
infrared K' band in two observing runs of 4 and 3 nights. 



respectively. The selection of these 24 observed sample 
galaxies depended solely on the allocation of telescope 
time; the galaxies cover the southern sky rather uniformly. 

By applying a V/Vmax completeness test (e.g., Davies 
1990; de Jong & van der Kruit 1994) we derived that the 
ESO-LV is statistically complete for diameter-limited sam- 
ples with D25 > l.'O. To check the completeness of our 
subsample, we calculated, based on a limiting diameter 
£>f5 > 2.'2, that I^/Kiax = 0.502 ± 0.253, which impHes 
statistical completeness. 

The near-infrared observations were obtained with the 
IRAC2B camera at the ESO/MPI 2.2m telescope of the 
European Southern Observatory (ESO) in Chile. The 
IRAC2B camera is equipped with a Rockwell 256x256 
pixel NICM0S3 HgCdTe array. For both observing runs, 
in July 1994 and January 1995, we used the IRAC2B 
camera with Objective C, corresponding to a pixel size 
of 0."491 (40 ^^m) and a field of view of 125" x 125". 

At both runs we used the K' filter available at 
ESO (central wavelength Ac = 2.15/im, bandpass AA = 
0.32/im). We chose to observe in K' rather than in if band 
(with Ac = 2.2/im, and AA = 0.40/im), since the K' band 
is almost as little affected by dust as the K band, but has 
a lower sky background (Wainscoat & Cowie 1992). 

We took sky images and object frames alternately, 
both with equal integration times (in sequences of 12 x 
10s), and spatially separated by ^ 5'. 

Supplementary observations in the Thuan & Gunn 
(1976) / band were obtained during a number of observ- 
ing runs at ESO. Most of the /-band observations were 
obtained with the Danish 1.54m telescope, equipped with 
a 1081x1040 pixel TEK CCD with a pixel size of 24^m 
(0.36 "/pix). The field of view thus obtained is 6'. 5 x 6.'2. 
The TEK CCD was used in slow read-out mode in order 
to decrease the pixel-to-pixel noise. The Thuan & Gunn 
(1976) /-band characteristics match those of a Johnson 
/filter (Buser 1978). For all observed galaxies we deter- 
mined the colour terms required for the calibration to the 
Cousins system using standard stars. 

Gaps in the observed sample were filled in by service 
observations with the Dutch 0.92m telescope, equipped 
with a 512x512 pixel TEK CCD. It has a pixel size of 
27 ^m (0.44 "/pix), corresponding to a field of view of 
3.'9 X 3.'9. 

Both telescopes were used in direct imaging mode, at 
prime focus. Details of the speciflc observations can be 
found in Table 0. 

During the reduction of the near-infrared observations, 
each sky frame was compared with the two sky frames 
taken nearest in time in order to detect stars in the sky 
frames. These stars were flltered out by using a median 
filter and thus the resulting cleaned sky images are very 
similar to the actual sky contributions. 

To circumvent the effects of bad pixels and to obtain 
accurate flatfielding we moved the object across the array 
between subsequent exposures. Therefore, for most galax- 
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Table 1. Log of the / and i^'-band observations 
Columns: (1) Galaxy name (ESO-LV); (2) Telescope used (Dan 
1.5 = Danish 1.54m; Dut 0.9 = Dutch 0.92m; ESQ 2.2 = 
ESO/MPI 2.2m); (3) Date of observation (ddmmyy) (4) Pass- 
band observed in; (5) Exposure time in seconds; (in the K' 
band the same integration time was spent to observe sky im- 
ages) (6) Seeing FWHM in arcsec 



Table 1. (Continued) 



Galaxy 


Telescope 


Date 


Band 


Exp. time 


Seeing 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


026G-06 


Dan 1.5 


090794 


/ 


2x900 


1.6 




ESQ 2.2 


110794 


K' 


2x12x10 


1.5 


033G-22 


Dan 1.5 


120194 


I 


2x900 


0.9 


041G-09 


Dan 1.5 


100794 


I 


2x900 


1.4 




ESO 2.2 


090794 


K' 


3x12x10 


1.5 


074G-15 


Dan 1.5 


110794 


I 


2x900 


1.3 


138G-14 


Dan 1.5 


080794 


I 


2x900 


1.6 


141G-27 


Dan 1.5 


090794 


I 


2x900 


1.2 




ESO 2.2 


090794 


K' 


2x12x10 


1.4 


142G-24 


Dan 1.5 


110794 


I 


2x900 


1.4 




ESO 2.2 


090794 


K' 


2x12x10 


1.4 


157G-18 


Dan 1.5 


100194 


I 


2x900 


1.2 




ESO 2.2 


240295 


K' 


2x12x10 


1.1 


201G-22 


Dan 1.5 


090194 


I 


2x900 


1.2 




ESO 2.2 


240295 


K' 


2x12x10 


1.1 


202G-35 


Dut 0.9 


061094 


I 


2x900 


1.2 


235G-53 


Dan 1.5 


120794 


I 


600 


1.5 




Dut 0.9 


170396 


I 


2x1200 


1.3 


240G-11 


Dut 0.9 


051094 


I 


2x1200 


1.6 






061094 


I 


2x1200 


1.6 


263G-15 


Dan 1.5 


090194 


I 


900 


1.0 






100194 


I 


600 


1.1 






120194 


I 


900 


1.0 




ESO 2.2 


240295 


K' 


12x10 


1.4 






250295 


K' 


12x10 


1.4 


263G-18 


Dut 0.9 


190194 


I 


2x900 


1.2 


269G-15 


Dut 0.9 


130396 


I 


2x1200 


1.3 


280G-13 


Dut 0.9 


140396 


I 


2x1200 


1.4 


286G-18 


Dan 1.5 


090794 


I 


2x900 


1.3 




ESO 2.2 


090794 


K' 


12x10 


1.6 






110794 


K' 


12x10 


1.2 


288G-25 


Dut 0.9 


051094 


I 


2x1200 


1.3 


311G-12 


Dan 1.5 


120194 


I 


2x900 


1.0 




ESO 2.2 


240295 


K' 


2x12x10 


0.9 


315G-20 


Dan 1.5 


090194 


I 


2x900 


1.7 




ESO 2.2 


240295 


K' 


2x12x10 


1.1 


321G-10 


Dut 0.9 


280493 


I 


300 


1.2 


322G-73 


Dut 0.9 


230396 


I 


2x1200 


1.2 


322G-87 


Dut 0.9 


270493 


I 


300 


1.3 


340G-08 


Dan 1.5 


100794 


I 


2x900 


1.3 


340G-09 


Dut 0.9 


140396 


I 


2x1200 


1.4 




ESO 2.2 


090794 


K' 


2x12x10 


1.5 


358G-26 


Dan 1.5 


120194 


I 


2x900 


1.0 


358G-29 


Dan 1.5 


100194 


I 


2x900 


1.3 




ESO 2.2 


260295 


K' 


2x12x10 


1.5 


377G-07 


Dut 0.9 


100396 


I 


2x1200 


1.2 


383G-05 


Dan 1.5 


080794 


I 


2x900 


1.3 




ESO 2.2 


090794 


K' 


2x12x10 


1.1 


416G-25 


Dan 1.5 


130194 


I 


2x900 


1.0 




ESO 2.2 


260295 


K' 


2x12x10 


1.4 


435G-14 


Dut 0.9 


201293 


I 


2x900 


1.0 




ESO 2.2 


240295 


K' 


2x12x10 


1.5 



Galaxy 


Telescope 


Date 


Band 


Exp.time 


Seeing 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


435G- 


•25 


Dut 0.9 


230493 


/ 


600 


1.7 








231293 


/ 


2x900 


1.3 








060194 


/ 


2x900 


1.2 






ESO 2.2 


240295 


K' 


2x12x10 


1.6 


435G- 


•50 


Dan 1.5 


100194 


I 


2x900 


1.1 


437G- 


•62 


Dut 0.9 


040195 


I 


2x1200 


1.0 






ESO 2.2 


260295 


K' 


2x12x10 


1.1 


444G- 


•21 


Dut 0.9 


230396 


I 


2x1200 


1.2 


446G- 


•18 


Dan 1.5 


090794 


I 


2x900 


1.1 






ESO 2.2 


090794 


K' 


2x12x10 


1.3 








120794 


K' 


2x12x10 


1.0 


446G- 


■44 


Dut 0.9 


140396 


I 


2x1200 


1.2 






ESO 2.2 


090794 


K' 


2x12x10 


1.4 


460G- 


•31 


Dan 1.5 


080794 


I 


2x900 


1.5 






ESO 2.2 


090794 


K' 


2x12x10 


1.5 


487G- 


•02 


Dan 1.5 


090194 


I 


2x900 


1.1 






ESO 2.2 


240295 


K' 


2x12x10 


1.1 


500G- 


■24 


Dan 1.5 


130194 


I 


2x900 


1.0 






ESO 2.2 


250295 


K' 


2x4x30 


1.6 


505G- 


■03 


Dut 0.9 


010496 


I 


2x1200 


1.3 


506G- 


■02 


Dut 0.9 


280493 


I 


300 


1.0 


509G- 


■19 


Dan 1.5 


100794 


I 


2x900 


1.0 






ESO 2.2 


100794 


K' 


2x12x10 


1.2 


531G- 


■22 


Dut 0.9 


270994 


I 


2x1200 


1.5 


555G- 


■36 


Dan 1.5 


130194 


I 


2x900 


1.1 


564G- 


•27 


Dut 0.9 


211293 


I 


2x900 


1.4 








221293 


I 


2x900 


1.3 






ESO 2.2 


240295 


K' 


2x12x10 


0.9 


575G- 


•61 


Dut 0.9 


230396 


I 


2x1200 


1.2 



ics mosaicing of cither 4 or 8 image frames was required to 
obtain complete galaxy images. The mosaicing was done 
by using common stars in the frames to determine the ex- 
act spatial offsets. In the rare case that no common stars 
could be determined, we used the telescope offsets as our 
mosaicing offsets. The overlapping area was used to deter- 
mine the adjustment of sky levels needed by means of a 
least squares fit. 

Bad pixels and bad areas on the array were masked out 
and not considered during the entire reduction process. 
Only after mosaicing was finished, the areas that still did 
not contain any valid observations were interpolated by 

a 2-dimcnsional linear plane fit (see Pclcticr [1993] for a 
detailed description of the reduction method used) . 

The calibration of the near-IR observations was done 
by using the SAAO/ESO/ISO Faint Standard Stars 
(Carter & Meadows 1995). We used the corrections pub- 
lished by Wainscoat & Cowie (1992) to transform the K' 
measurements to the K band. The accuracy of the K'- 
band zero-point offsets we could reach was ~ 0.08 mag at 
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both observing runs. The hmiting factors here were flat- 
fielding errors. 

The /-band images were reduced foUowing standard re- 
duction procedures (see de Grijs & van der Kruit 1996); for 
the cafibration of these observations Landolt fields were 
used (Landolt 1992). The /-band calibration could be done 
to an accuracy of ~ 0.03 — 0.05 mag, depending on the 
telescope and observing run. 

Both our /-band observations and the /sT'-band data 
were taken at photometric (parts of) nights. 

2. 2. Vertical profiles 

We extracted vertical luminosity profiles at a number of 
positions along the major axes of the sample galaxies. 
A semi-logarithmic binning algorithm was applied to the 
galaxies both radially and vertically, in order to retain an 
approximately constant overall signal-to- noise (S /N) ratio 
in the resulting vertical profiles. Since the most significant 
differences between our models become clear at small z, 
no vertical binning was applied close to the galaxy planes. 

We rejected those profiles with low S/N ratios (gen- 
erally the outermost profiles) and those that were clearly 
affected by artifacts in the data or foreground stars. For 
all galaxies we have been able to sample the vertical light 
distribution at various positions along the major axis out- 
side the region where the bulge contribution dominates. 
For galaxies of types T > —1.0 we could determine these 
distributions for at least 4 of these independent positions. 

The positions of the galaxy planes were determined by 
folding the vertical profiles and under the assumption of 
symmetrical light distributions with respect to the planes, 
in the near- infrared /C'-band observations. 

In Fig. 0we present the /-/f colour maps of the galaxies 
discussed in this paper. 

2.3. Comparison to published luminosity profiles 

Surface brightness profiles of edge-on galaxies observed 
with modern detectors are scarcely available in the recent 
literature. 

In the near-infrared K band we compared our observa- 
tions of the large southern edge-on galaxy ESO 435G-25 
with those of Wainscoat et al. (1989), obtained using a 
raster scan technique with an aperture of 5". Although 
our observations of ESO 435G-25 are of a much higher 
quality and were taken with a much higher resolution, 
we find remarkably good agreement between Wainscoat et 
al.'s (1989) and our /f-band observations, as can be seen 
in Fig. ||. We have extracted vertical profiles from our cali- 
brated /T-band image at exactly the same positions and in 
the same way as was done by Wainscoat et al. (1989). The 
distinct drop in the difference profiles at the position of 
the galaxy plane is caused by the difference in resolution 
between both data sets. 
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Fig. 2. Comparison of K-hand vertical profiles obtained by 
Wainscoat et al. (1989) and our observations. To obtain these 
profiles we averaged strips of 7. "5 wide. 



A large and homogeneous set of /-band observations 
of southern Sb-Sc galaxies has recently been published by 
Mathewson et al. (1992) and Mathewson & Ford (1996). 

Azimuthally averaged luminosity profiles were ob- 
tained by fitting ellipses to the galaxy isophotes, whose in- 
tensity, ellipticity and position angle were allowed to vary 
with each ellipse. Bad pixels, cosmic rays and foreground 
stars were masked out, so that these would not affect 
the results from the ellipse-fitting routine. Although this 
method works sufficiently well for low and moderately- 
inclined galaxies, when dealing with highly-inclined or 
edge-on galaxies the ellipse fitting is severely infiuenced by 
the presence of a central dust lane and the non-elliptical 
outer galaxy isophotes. Unfortunately, since Mathewson 
et al. (1922) and Mathewson & Ford (1996) did not tabu- 
late the ellipticities nor the position angles used for the 
individual ellipses obtained for each galaxy, we can at 
best compare azimuthally averaged profiles which were 
obtained with the same free parameters. A comparison 
between the azimuthally averaged /-band luminosity pro- 
files of Mathewson et al. (1992) and Mathewson & Ford 
(1996) and those obtained from our observations is shown 
in Fig. I 

In general, we find that the differences between our and 
Mathewson's measurements are small, although clear de- 
viations are appreciated in a number of cases. In particular 
for those galaxies for which the difference between our and 
Mathewson's profiles is relatively large (e.g., ESO 138G- 
14), we used observations obtained on different nights or 
with a different telescope to check our results. It was found 
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Fig. 1. Calibrated I-K colour maps of our sample galaxies, displayed on the same scale. Each colour map has been inserted 
in a rectangular frame of size 80" x 200". Galaxy images were rotated from their original sky orientations so that the major 
axes lie vertically. The corresponding calibrated /-A" colours are indicated by the gray-scale bar at the bottom of the figure, (a) 
ESO026G-06; (h) ESO041G-09; (c) ES0141G-27; (d) ES0142G-24; (e) ES0157G-18; (f) ESO201G-22; (g) ES0263G-15; (h) 
ES0286G-18; (i) ES0311G-12; (]) ESO315G-20; (k) ESO340G-09; (I) ES0358G-29; (m) ESO383G-05; (n) ES0416G-25; (o) 
ES0435G-14; (p) ES0435G-25; (q) ES0437G-62; (r) ES0446G-18; (s) ES0446G-44; (t) ESO460G-31; (u) ESO487G-02; (v) 
ESO500G-24; (w) ESO509G-19; (x) ES0564G-27. 



that the features shown in Fig. ^ can be reproduced to 
within the observational errors. The main cause of devia- 
tions between Mathewson's and our profiles, in particular 
at small semi-major axis radii, is the unpredictable influ- 
ence of dust, which greatly affects the ellipse fitting in the 
inner galaxy regions. 

From this comparison to previously published surface 
brightness profiles, we conclude that our observations re- 
produce both the photometric zero points and the be- 
haviour of the galaxy light as a function of position across 
the galaxy disk to within the observational errors. 

2.4. Extinction correction 

As was shown by Aoki et al. (1991), the inner contours 
of the if-band image of NGC 891 are asymmetric. This 
shows that absorption is not completely negligible, even 
at these near-infrared wavelengths. 

To correct for this small absorption effect at K, we used 
the J-Ti' colour index as an extinction indicator, following 
Knapen et al. (1995). We assumed that deviations from 
an average I-K colour are due to dust only and that the 
Galactic extinction law applies also in external galaxies 
(see, e.g., Jansen et al. 1994). Although the exact use of 
this law requires a detailed knowledge of the geometry of 
the mixture of dust and stars, any errors caused by our 
assumptions are of second order, since the K extinction 
due to dust is small. 

As an example, in NGC 891, Aoki et al. (1991) pointed 
out that the southwest part of the if-band image is rather 
patchy along the major axis, compared to the northeast 
part. The most likely explanation for this patchiness is 
that it is due to dust associated with spiral arms, as can 
be deduced from colour images. 

In our assumptions, we ignore two kinds of system- 
atic errors, as was pointed out by Knapen et al. (1995): a 
contribution of smoothly distributed dust, which does not 
alter the galaxy's morphology, and effects due to popula- 
tion changes, which are believed to be small. In general I-K 
colours of stellar populations are very similar. For exam- 
ple stellar population models (e.g., Vazdekis et al. 1996) 
only show a small range for models of different ages and 
metallicities. Although population gradients probably are 
present, the errors made by not taking them into account 
are likely to be so small that it is better to apply this 
extinction correction than not to apply it. 



Summarized, we corrected our profiles as follows: 

i^corr = iCobs - ({I - KUs -{I- i^)avg)) , (5) 

where -ftTcorr and iCobs are the corrected and observed pro- 
files, (/-ii:)obs and {I-K) avg are the observed and mean 
{I~K) colours, and Aj and Ak are the extinction in / and 
K, respectively. We used Aj /Ak = 4.30 (Ricke & Lebofsky 
1985). This extinction correction amounts to 0.55 mag in 
K at maximum, thus showing that we are dealing with op- 
tically thin regions in our galaxies. Therefore, by applying 
this correction, the errors that are caused by the incorrect 
underlying assumption that the dust is distributed in a 
foreground screen, are small. 

2. 5. A generalized family of density laws 

By analyzing these K'-hand images of our sample of edge- 
on disk galaxies we should be able to distinguish statisti- 
cally between the various models for the vertical luminos- 
ity and mass density distribution. 

As it seems reasonable to take the isothermal and ex- 
ponential distributions as the two extremes, van der Kruit 
(1988) proposed to use the family of density laws 

K{z) = 2-^/"Ko sech2/"(nz/2zo), (n > 0) (6) 

where K(z) is the observed ii'-band vertical density pro- 
file, Kq is the extrapolated outer surface density in the 
galaxy plane, z is the distance from the plane, and zq is 
the vertical scale parameter. The isothermal model is the 
extreme for n — 1, and the exponential is the other ex- 
treme for n = oo. For comparison, for the isothermal sheet 
Zq = 2hz, where hz is the exponential vertical scale height. 
Therefore, the sharpness of the peak in a luminosity pro- 
file is determined by the exponent 2/n of this family of 
density laws (^. 

In Fig. ^we plot a few model luminosity distributions, 
with the exponential and the isothermal functions as the 
two extremes. We have adopted identical central surface 
brightnesses and vertical scale heights for each of these 
models. 

3. Results 

3.1. Near-infrared profiles 

In order to quantitatively distinguish between the various 
models, we tried to fit the extinction-corrected vertical 
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Fig. 3. Comparison between azimuthally averaged /-band profiles published by Mathewson et al. (1992) and Mathewson & 
Ford (1996), and those obtained from fitting ellipses to our observations. 



8 



R. de Grijs et al.: The 2-structure of disk galaxies 




z height (z„) 



Fig. 4. 

The family of density (luminosity) laws ^ with the isothermal 
(2/n = 2.0) and the exponential (2/ra — 0.0) distributions as 
the two extremes. The difference in 2/n between two succes- 
sive model distributions is O.f 25; for clarity, the sech(2) model 
distribution (2/n = 1.0) is shown as the dashed profile {fio is 
the central surface brightness). 



light profiles in the K band with the extrapolated central 
surface brightness, the vertical scale parameter and the 
exponent 2/n as free parameters. To check the validity 
of the best-fitting 2/n thus obtained, we also tried to fit 
the vertical luminosity profiles with a continuous series 
of fixed values for the 2/n parameter. We compared the 
minimum in the distribution as a function of 2/n with 
the best-fitting n value resulting from the three-parameter 
fits. From this comparison, we concluded that resulting 
2/n values from the three-parameter fits could reproduce 
the minima in the distribution with sufficient accuracy 
to be reliable indicators of the profiles' shapes. 

To show the accuracy of our fits, in Fig. |^ we have com- 
piled all our vertical i^'-band profiles. We show the best- 
fitting model luminosity distributions obtained with the 
three-parameter fitting routine for those profiles that are 
sufficiently symmetric. The vertical luminosity profiles for 
each galaxy are shown in two separate panels, represent- 
ing either side of the galaxy with respect to the galaxy 
center; for comparison, the central profiles are shown in 
each panel. In Table || we have listed the positions along 
the major axis at which the vertical luminosity profiles 
were extracted. 



Table 2. Positions along the major axis at which the 
vertical profiles of Fig. 5 were extracted. 

Columns: (1) Profile number (on either side of the galaxy cen- 
ter; 1 — central profile) ; (2) center position of radially averaged 
bin along the major axis (arcseconds) ; (3) bin width (arcsec- 
onds) . 



Profile number 


Center position 


Bin width 


(1) 


(2) 


(3) 


1 


0.0 


10.8 


2 


11.8 


12.9 


3 


26.0 


15.4 


4 


43.1 


18.7 


5 


63.5 


22.3 


6 


88.0 


26.6 



In Fig. |6|we present the distribution of the best-fitting 
values for the exponent of the family of luminosity laws 
(||) for our total sample of K'-haml images. 

We notice that, apart from the rather unpredictable 
behaviour in the central areas, where the disk luminosity 
is contaminated by a non-negligible bulge light contribu- 
tion, the vertical luminosity profiles generally exhibit little 
variation of 2/n with position along the major axis. As has 
been shown by Andredakis et al. (1995), the early-type 
bulges are characterized by very steeply peaked luminos- 
ity profiles, whereas in the later-type galaxies the bulge 
profiles become closer to exponential. Therefore, we ar- 
gue that the dominant bulge light contribution causes an 
increase of the cuspiness of the profiles, especially in the 
earlier- type galaxies. 

In all panels of Fig. |^, we have indicated an estimated 
mean level of the best- fitting 2/n values by a dash-dotted 
line. These values are given in Table ^. In determining 
these mean levels, we have excluded both those fits that 
were substantially affected by bulge light contamination 
and those fits for which the resulting values of the 
minimization routine were excessively high, due to either a 
low signal-to-noise ratio or foreground stars superimposed 
on the galaxy. Especially in the (radially) outer parts of 
the galaxy disks our K' observations are rather noisy and 
therefore unreliable. 

In general, we find that the mean levels for the sharp- 
ness of the K' band luminosity peaks indicate that the 
vertical luminosity distributions are more peaked than ex- 
pected for the intermediate sech(z) function, proposed by 
van der Kruit (1988), though less peaked than exponen- 
tial. 

To study the variation of the cuspiness of the vertical 
profiles as a function of position along the major axis, we 
have averaged the data points of our A''-band observations 
in radial bins of 0.5 scale length, both for our total sample 
and for subsamples in certain ranges of revised Hubble 
type, T. The results are shown in Fig. |^. The central values 
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Fig. 6. 

Best-fitting exponents determining the sharpness of the peak in the galaxies' vertical luminosity profiles. The galaxies are 
ordered by revised Hubble type, which is shown in the upper right corner of each panel. The dotted lines indicate the levels 
corresponding to an exponential profile (2/n = 0), an isothermal profile (2/n = 2), and the intermediate soch(«) distribution (2/n 
= 1). The dash-dotted line represents the mean best-fitting value for the galaxies' disks, obtained by averaging the filled data 
points; the open circles were not taken into account when determining these levels, because of either bulge light contamination 
or an excessively high value. The positions along the major axis at which the luminosity profiles were extracted are expressed 
in units of their if-band radial scale length, hr,K- 
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Table 3. Best-fitting 2/n values 

Columns: (1) Galaxy name (ESO-LV); (2) Revised Hubble 
Type (T); (3) and (4) Mean 2/n values and the (1 a) error 
obtained from the X'-band observations. 



Galaxy 


Type 




± 


(1) 


(2) 


(3) 


(4) 



0.10 
0.17 
0.10 
0.10 
0.09 
0.24 
0.13 
0.30 
0.15 
0.34 
0.16 
0.08 
0.04 
0.21 
0.27 
0.12 
0.16 
0.08 
0.14 
0.21 
0.13 
0.13 
0.22 
0.16 



for 2/n, (2/n)o, were used for normalisation; the positions 
along the major axis are expressed in units of the iiT-band 
radial scale lengths, hr,K- We determined the scale lengths 
by fitting ellipses to the two-dimensional galaxy isophotcs 
in the regions away from the central dust lane, roughly 
between 1 and 4 radial scale lengths, depending on the 
bulge influence. 

We chose to look at the galaxies in the type range 
—2.0 < T < 1.0 (4 galaxies) because of their clear bulge 
contribution, and those in the range 4.0 < T < 7.0 (13 
galaxies), since they are well-behaved late- type galaxies 
without prominent bulges. The intermediate type range, 
with 1.0 < T < 4.0, is shown as well. For these galaxies the 
bulge influence is generally small, though not negligible. 
In order to reduce the noise in the earliest-types bin, the 
radial binning was done in intervals of 1.0 scale length. 

We notice that the distribution of the sharpness of the 
vertical profiles remains, within the errors, constant as 
a function of position along the major axis, irrespective 
of galaxy type. Some of the (radially) outermost profiles 
show a slight increase of 2/n, corresponding to a rounder 
vertical profile, but in those regions the number of useful 
profiles is small and hence the errors are large. The inner- 
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most profiles, especially those in the earliest-types bin, are 
affected by bulge light and should therefore be left out of 
the analysis. 
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Fig. 7. 

Averaged distributions of the sharpness of the vertical profiles 
as a function of position along the major axis, both for the 
total /f'-band sample and for specific bins in revised Hubble 
type, T. 



3.2. Effects of seeing 

Since we are particularly interested in the behaviour of 
the vertical light distributions in edge-on disk galaxies at 
small z, we have to be careful to account for the effects 
of atmospheric smearing. No matter what the exact light 
distribution is, atmospheric seeing causes the profiles to 
show a flat-topped distribution in the galaxy planes. In 
order to circumvent this problem, wc have to convolve our 
models with the seeing proflle. Assuming that the actual 
point spread function can be adequately described by a 
Gaussian, the seeing convolved model profile is: 

/■OO 

K^{z) = a-^e-yc^{-z'^/2a^) / K{x)lQ{xz/a^) 

Jo 

X exp(-a;^/2CT^)a;da;, (7) 

where Kc{z) and K(z) are the corrected and intrinsic 
model intensity profiles, a is the dispersion of the seeing 
Gaussian and Iq is the zeroth-order modified Bessel func- 
tion of the first kind (Pritchet & Kline 1981; Andredakis 
& Sanders 1994). 
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The importance of atmospheric smearing on a partic- 
ular vertical profile depends on the ratio of the seeing 
FWHM and the vertical scale height of the profile. The 
effects are most noticeable for (nearly) exponential pro- 
files. In the case of an exponential light distribution, the 
seeing convolution becomes non-negligible for those ob- 
servations with seeing FWHM > 0.6 h^, see Fig. ^a. In 
the present study of near-infrared observations, the see- 
ing FWHM varies between 0.12 and 0.57 h^. Therefore 
we conclude that the (variation of the) cuspiness of the 
observed vertical profiles is real rather than an artifact of 
the method applied. This conclusion is supported by the 
fact that, in some of the earliest-type galaxies, we see a 
variation of the 2/n parameter along the major axis due 
to bulge light contamination, causing more sharply peaked 
vertical profiles, which is likely to be real. 



-Taj 



- 3ech(z) 
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Fig. 8. 

(a) The influence of atmospheric smearing on an exponential 
vertical luminosity proflle. The difference between the true and 
the observed distributions is shown as a function of the seeing 
FWHM, in units of the profile's scale height; (b) Distribution 
of absolute magnitude versus 2/n, the typical error size is in- 
dicated in the lower left corner. 



3.3. Comparison with I-band observations 

To be able to compare our results with those published 
previously, we applied this method to our total sample 
of /-band observations, thereby using the same fitting 
ranges as for the i^'-band data. It is immediately obvi- 
ous that the 7-band results are much noisier than those 
obtained in the K' band, largely due to the presence of 
clear dust lanes along the galaxies' major axes. In those 
galaxies that show a well-behaved dependence of the 2/n 
value as a function of position along the galaxy's ma- 
jor axis, we cannot distinguish statistically between the 
best-fitting 2/n values in the / band and those obtained 
in the K' band, within the errors. This is made clear 
in Fig. H, where we compare the distributions of best- 
fitting 2/n levels between the two passbands. If we ap- 
proximate the distributions of 2/n in the / and K' bands 



by a Gaussian, we find {2/n)K' — 0.538, ctk' = 0.198 and 
(2/n)/ = 0.528, cr/ = 0.202. However, K'-hand observa- 
tions are clearly preferable to determine the sharpness of 
the peak of the vertical luminosity profiles unambiguously, 
because of their relative insensitivity to contamination by 
galactic dust. In the presence of a central dust lane, i.e. 
a deep trough in the vertical light distribution, one has 
to remove the central regions from the fit to get mean- 
ingful answers. While doing this, one often ends up only 
with the exponential outer parts of the vertical distribu- 
tion, forcing the fitting routine to yield a more sharply 
peaked solution. Although we tried to fit the vertical pro- 
files over the same range in z for both the / and K'-hand 
observations, this was not always possible due to the dis- 
turbing effects of dust in the /-band data. Therefore the 
best-fitting /-band 2/n values are not very reliable. In K' 
the best-fitting 2/n value is hardly affected by a different 
fitting range, since the fit is dominated by the inner data 
points, which represent the luminosity distribution in the 
galaxy plane. 

I and K'-band data available 
Only /-band data available 





0.5 1 1.5 



Fig. 9. 

Comparison between the best-fitting exponents (2/n) obtained 
from both K'-hand observations (left panel) and /-band obser- 
vations (right panel). In the right panel, we show the subsample 
corresponding to the K' subsample shaded differently. Statis- 
tically, these two distributions are indistinguishable. 



From Fig. ^ it can be seen that the median of the best- 
fitting 2/n value lies in between the exponential (2/n = 
0.0) and the intermediate sech(z) model (2/n — 1.0). It 
is therefore not surprising to note that in previously pub- 
lished papers by, e.g., Wainscoat et al. (1989) for the large 
southern edge-on galaxy IC 2531, and Gilmore & Reid 
(1983) for the Galaxy, the exponential and the sech(2;) 
models were found to fit the data equally well (van der 
Kruit 1988). 

When studying the relation between the sharpness of 
the peaks and various global parameters for our sample 
galaxies, little correlation was detected, although it seems 
that the smallest galaxies generally show a more sharply 
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peaked vertical distribution than the larger galaxies. In 
other words, for those galaxies with the faintest absolute 
magnitudes (M/ > —20 or Mk > —22), we find a lack 
of "rounder" profiles (2/n > 0.6) compared with those 
galaxies with brighter absolute magnitudes, as can be seen 
in Fig. 

Although we do not find clear correlations with global 
galaxy parameters, we argue that the radial distribution 
of 2/n is universal (i.e., we do not find any significant 
variations as a function of galaxy type), exhibiting only 
local variations due to, e.g., residual dust contamination 
(see, e.g., Wainscoat et al. 1989; Aoki et al. 1991). 

3.4^. Extinction analysis 

In Fig. |l^ we show the central profiles of all sample galax- 
ies in both / and K' band, and their corresponding I-K 
colour, to give an indication of the importance of contam- 
ination by dust in the /-band observations. 

The vertical colour profiles seem to indicate strong ex- 
tinction in the plane of the galaxy, and no extinction in 
the outer parts. We therefore define the constant colour in 
the outer parts to be the colour of the stars in the galaxy, 
and the maximum I-K difference to be the I-K colour 
excess, E{I — K). 

In Table ^ we give, on the first line, crude estimates for 
the maximum (edge-on) /-if colour excess, E{I — K)■a^a:^^ 
in the galaxy planes, at various positions along the galax- 
ies' major axes. We obtained these estimates by averaging 
our vertical colour profiles on both sides of the galaxy 
centers. In deriving these colour excesses we assume that 
there is little or no dust mixed in with the stellar popu- 
lation outside the dust lane region, as seems likely from 
the smooth vertical colour profiles. From a comparison 
with published colours of moderately inclined Sc galax- 
ies, Kuchinski & Terndrup (1996) have shown that there 
is indeed little or no reddening away from the dust lane. 
Secondly, we have neglected any change in intrinsic colour 
due to a stellar population gradient by assuming that the 
intrinsic colours of the dominant stellar population in the 
dust lane region are identical to those away from the dust 
lane. Kuchinski & Terndrup (1996) argue that a popula- 
tion gradient would result in a very small colour change 
compared to the total colour change across the dust lane 
and that we can therefore assume that the colour exces is 
largely due to dust extinction. 

By using the Galactic extinction law we derive that 
E{I — K) = 0.77/1/. 900, the maximum (edge-on) /-band 
extinction. From this, one can derive the optical depth, 
and the optial depth in V: Ty_9o° ~ 2t/_9oo . 

In all our sample galaxies, E{I — K)-^g_^ remains con- 
stant or decreases as a function of projected galactocentric 
distance, which could be a characteristic of an exponen- 
tially decreasing dust contribution. Jansen et al. (1994) 
find, for a small sample of highly inclined spiral galaxies. 



that the maximim extinction in the dust lane (^v.max) 
decreases rapidly with increasing galactocentric distance. 

For each position for which we give the maximum /- 
K colour excess, we also tabulate the inferred maximum 
/-band optical depths, r/^max, assuming a uniform mix- 
ture of dust and stars in the galaxy planes, taken from 
Walterbos & Kennicutt (1988): 

A, = -2.51ogio(^^) (8) 

For the same optical depth, the uniform mixture of 
dust and stars causes less /-band extinction, A/, than the 
classical foreground screen model, because part of the ex- 
tinction lies behind the source. In Table ^ we have also 
listed the inferred face-on optical depth estimates. For 
face-on galaxies we have assumed the galaxies to be 9 
times less opaque than in the highly-inclined case, which 
reflects the difference in path length through the disk at 
both positions. Of course, the face-on optical depth can 
only be estimated in a rough way; in reality the opti- 
cal depth depends on the geometry of the dust and stars 
(which we have assumed to be uniformly mixed) , the filling 
factor of the disk components (in our approximation we 
have neglected the possibility of a patchy dust distribution 
and the presence of spiral arms), and the major-to- minor 
axis ratio, which does not need to be 9 for all our sample 
galaxies. 

In Fig. |l^a we present the mean values for the maxi- 
mum /-band edge-on optical depth as a function of Hubble 
type and for different positions. The scatter in this figure 
is large due to the small number of data points available, 
but also because of the uncertainties in the galaxy clas- 
sification inherent to edge-on galaxies. Therefore, we can 
only draw qualitative inferences from the behaviour of the 
optical depth as a function of galaxy type. This behaviour 
suggests an increasingly important dust contribution from 
the lenticular and early spiral galaxies towards later types, 
although for the latest galaxy types the dust content seems 
to diminish relative to the intermediate types. 

This result is in qualitative agreement with that of 
Peletier et al. (1994), and Peletier & Balcells (1997) (see 
also Peletier & Balcells 1996), who studied the type de- 
pendence of the R and / vs. K' radial scale length ratios 
of a sample of some 70 disk galaxies, both face-on and 
edge-on. From Fig. |ll|& it is clear that the radial colour 
gradients, indicated by deviations from unity of this ratio, 
are smallest for the early-type sample galaxies, whereas in 
the later types they vary considerably. 

Peletier et al. (1994) have shown that scale length ra- 
tios due to stellar population changes are of order 1.1-1.2 
in the blue - near-infrared range; in the / vs. K' range 
this contribution is likely to be less. Therefore, the ob- 
served scale length ratios largely represent the galaxies' 
dust content. 

Although the scale length ratios indicate a more or less 
constant dust content for galaxy types later than about 
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Table 4. /-band extinction estimates in the galaxy 
planes 

Columns: (1) Galaxy name; (2)-(4) first line: Extinction es- 
timates at the galaxy center, at 1 and at 2 radial (A''-band) 
scale lengths. The accuracy of the numbers is indicative of the 
quality of the colour profile; second line: Inferred maximum 
/-band optical depths, r/,max, assuming a uniform mixture of 
dust and stars (edge-on and face-on values). 



T = 3 (Sb), the data suggests that the ratios decrease 
towards later types, as is also found from the optical depth 
measurements. 
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Fig. 11. 

(a) Maximum /-band optical depths as a function of galaxy 
type. The errors are statistical errors; where no error bars are 
given, only one data point was available in the specific type 
bin, indicated by the bars; (b) Type dependence of the / vs. K' 
scale length ratio. The filled circles are from the observations 
presented in this paper, the open circles from Peletier et al. 
(1994), and the crosses from Peletier & Balcells (1997). 



4. Discussion 

Main observational results 

As shown in the previous section, our main observational 
results are the following: 

— The mean levels for the sharpness of the iiT'-band lu- 
minosity peaks indicate that the vertical luminosity 
distributions are more peaked than expected for the 
intermediate sech(2) distribution, but rounder than ex- 
ponential. 

— In the majority of our sample of edge-on disk galaxies, 
we find that the sharpness of the peak of the vertical 
profiles, characterized by the exponent 2/n of the gen- 
eralized family of fitting functions (^) (van der Kruit 
1988), varies little with position along the major axis. 
This result is independent of galaxy type. 

— Due to the unpredictable contamination by in-plane 
dust, /iT'-band observations are preferred to /-band 
measurements, although statistically we cannot distin- 
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Fig. 10. 

Central 7 and if-band vertical profiles for the total sample, and their corresponding I-K colours 
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guish between the results obtained in cither of the two 
passbands. In general, the /-band determinations of 
the best-fitting Ij n values are much noisier than those 
obtained in the near-infrared. 

The behaviour of the maximum /-band optical depth 
as a function of galaxy type suggests an increasingly 
important dust contribution from the lenticular and 
early spiral galaxies towards later types, although for 
the latest galaxy types the dust content seems to di- 
minish relative to the intermediate types. 
For those galaxies with the faintest absolute magni- 
tudes, we find a lack of "rounder" profiles compared 
to those galaxies with brighter absolute magnitudes. 



^.2. Projection effects 

In interpreting the observed surface brightness distribu- 
tions as surface density representations, we should be cau- 
tious: different stellar populations may have both differ- 
ent mass-to-light ratios and different velocity dispersions. 
Moreover, the observed galaxy surface brightness distri- 
butions are line-of-sight integrations of the light contribu- 
tions and therefore include contributions from very differ- 
ent locations in the disk. However, as noted by, e.g.. Dove 
& Thronson (1993), these effects are normally considered 
to be small. 

Kylafis & Bahcall (1987) showed that the main effects 
of scattering are the reduction of the extinction in the dust 
lane due to forward scattering and an apparent thickening 
of the central layer of high luminosity. This would not lead 
to a more sharply peaked profile if the underlying light (or 
density) distribution were isothermal. Moreover, scatter- 
ing is not expected to be important in the near-infrared 
for edge-on galaxies (Kuchinski & Terndrup 1996). 

Since the differences between the models discussed 
here are small, we have investigated the effects of pro- 
jection. Projection effects lead to an apparent thickening 
of the high-luminosity midplane of the galaxy. Deviations 
from an inclination of 90° cause a slight increase of the 
scale height, as can be seen in Fig. |l^. The effect is com- 
parable for the exponential and the isothermal models. 
We are confident that none of our sample galaxies have 
inclinations lower than 87°. 

The effect on the 2/n values of deviations from a 90° 
inclination is shown in Fig. ^|b for an exponential dis- 
tribution at the galaxy center, at 1 and at 2 radial scale 
lengths, respectively. For exponential distributions the ef- 
fects of projection are most noticeable, compared to more 
flattened vertical luminosity distributions. In fact, the ef- 
fects of projection can be so large, that all our galaxies 
can have intrinsically exponential vertical surface bright- 
ness distributions. 
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Fig. 12. 

Importance of projection effects on the results obtained. An in- 
trinsic axis ratio of (^) — 0.11 was assumed (Guthrie, 1992). 
Projection effects are siiown for a vertical surface brightness 
profile at the galaxy center (short-dashed line) , at 1 radial scale 
length (dotted line) and at 2 radial scale lengths (long-dashed 
line), (a) Projection effects cause an increasing scale height 
with increasing deviations from a 90° inclination; (b) The in- 
fluence of projection effects on the 2/n parameter. 



4-. 3. Disk heating 

In this paper we show that the universal vertical structure 
that is observed away from the galaxy planes is maintained 
even when going all the way down to very low z distances. 
Models that aim to explain this observational fact should 
also be able to account for intrinsically exponential verti- 
cal surface brightness or density distributions. 

As Dove & Thronson (1993) warn, if the vertical distri- 
bution of stars cannot be fit well by the isothermal sheet 
approximation, which model is based on physical princi- 
ples, then rather than invoking an arbitrary alternative 
function, a more physical function should be that of a 
nonisothermal stellar distribution. A nonisothermal distri- 
bution of stars can be represented by a linear combination 
of isothermal components (Oort 1932). 

It is possible that once a coeval population of stars 
has been formed, they do not (or only over a long period 
of time) interact with other components and are there- 
fore quasi-independent isothermal components (Dove & 
Thronson 1993). Such a model would be physically re- 
alistic, since if stars have formed at different times with 
different velocity dispersions and have not yet reached an 
equilibrium state, or if the stars are dynamically heated 
during their evolution, the resulting distribution of stars 
would not be well approximated by a single isothermal 
model. From observational evidence we know that stel- 
lar subpopulations of different ages have different velocity 
dispersions. Star formation is generally believed to be a 
continuous process. Therefore, a more accurate model for 
a galaxy disk is a superposition of a very large number of 
components, or, as Kuijken (1991) proposed, an integral 
representation. 
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The observations presented in this paper indicate that 
the physical process responsible for the vertical luminosity 
profiles is probably both global and universal in nature, 
since we do not find any significant radial variations of the 
cuspiness of the profiles, nor large variations as a function 
of galaxy type. 

This puts interesting constraints on the dominant ver- 
tical heating mechanism in galaxy disks. 

Although the dominant perturbations in a disk are the 
spiral waves, they have a close relationship with the giant 
molecular clouds (GMCs). Juhan & Toomre (1966) and 
Julian (1967) suggested that the GMCs must acquire large 
wakes of material, thus increasing the effective mass of the 
combined spiral and GMC perturbation. These wakes can 
be very strong, but will also spread over a large area of 
the disk, which makes it more difhcult to assess the impor- 
tance of such wakes in enhancing the scattering efficiency. 
Moreover, at present there is no satisfactory explanation 
as to how disk heating, the rate of which must vary greatly 
with radius from the observed distribution of GMCs in our 
own and other galaxies, can naturally lead to a global and 
universal vertical density distribution and a constant scale 
height with galactocentric distance (i.e., zq is independent 
of R) (e.g., Jenkins 1992). 

4-4- Isothermal versus exponential distributions 

Jenkins (1992) finds, that his model disk heating process, 
i.e., combined spiral and GMC perturbations enhanced by 
disk accretion, together with constant star formation, al- 
ways leads to a closely isothermal stellar population. Wie- 
len (1977) reached a similar conclusion based on observa- 
tional data. 

Burkert & Yoshii (1996) show, based on realistic hy- 
drodynamical calculations of disk evolution processes, 
that - if one starts from a non-equilibrium gaseous state 
- the final vertical stellar density profile depends strongly 
on the initial distribution of the protodisk gas, as opposed 
to the GMC heating process described in the previous sec- 
tion. 

On the other hand, if they assume that the gas set- 
tles into isothermal equilibrium prior to star formation 
and gas cooling, then always an exponential density pro- 
file is formed, although the vertical scale height increases 
as a function of decreasing surface brightness. In fact, in 
de Grijs & Peletier (1997) we presented the results of a 
detailed study of the vertical scale height distributions in 
the present sample, for which we found an increasing scale 
height with galactocentric distance, in particular for the 
earlier-type galaxies. An interesting result from the calcu- 
lations of Burkert & Yoshii (1996) is that when the ratio 
of the star formation time scale to the cooling time scale 
lies in the range between ^ 0.3 and 3, the vertical stellar 
density distribution becomes exponential, independent of 
the free parameters in their modeling and also indepen- 



dent of the initial (isothermal) disk temperature and the 
initial surface density. 

Therefore, the process of crucial importance is that the 
SFR is adjusted sooner or later to balance with the local 
cooling rate (Burkert & Yoshii 1996). 

Just et al. (1996) find that if the SFR decreases with 
time, exponential luminosity profiles also grow naturally. 
They state that since the mass-to-light ratio of a stellar 
population increases with age, an exponential luminosity 
profile corresponds to a density profile that is slightly flat- 
tened to the galaxy plane. However, an isothermal density 
distribution is too thick to explain exponential light pro- 
files. For a constant SFR an obvious luminosity excess 
near the plane would show up in the optical bands, when 
assuming a heating mechanism of the type observed in the 
solar neighbourhood. 

The fact that we observe a more strongly peaked verti- 
cal light distribution than a sech(z) model in all our sam- 
ple galaxies, independent of galaxy type, indicates that 
the process at work here is a process intrinsic to the disks 
themselves, rather than a type-dependent mechanism. The 
variations in the cuspiness of our profiles along the galax- 
ies' major axes are probably due to some local mechanism, 
e.g. the contamination by residual dust. Although we ob- 
serve a similar behaviour in SOs as in later-type galaxies, 
this does not necessarily mean that they all possess young 
populations, although the young population contributes 
also in the near-infrared. The dominant stellar population 
in K' is the old population, with ages of several Gyrs, 
which is likely present in both early and late-type galax- 
ies. Therefore our observational result of a universal verti- 
cal density profile in the near-infrared is not incompatible 
with no type dependence. 

Finally, the fact that we observe a lack of "rounder" 
profiles in the smaller galaxies compared with the larger 
ones, may indicate that we are hindered by an underly- 
ing dust component, which is concentrated towards the 
galaxy planes, and more extended in the larger galaxies 
than in the smaller ones, at least to an outside observer. 
It may be that in the smaller galaxies this dust component 
affects relatively fewer data points than in the larger ones, 
thus causing a bias towards more sharply peaked vertical 
profiles. 

5. Summary and Conclusions 

In this paper we have studied the sharpness of the peaks 
of the vertical luminosity and density distribution in a 
statistically complete sample of edge-on disk galaxies. The 
results obtained in this paper are based on near-infrared 
if'-band observations. The main results obtained in this 
study are the following: 

— In the majority of our sample of edge-on disk galaxies, 
we find that the sharpness of the peak of the vertical 
profiles, characterized by the exponent 2/n of the gen- 
eralized family of fitting functions (^) (van dcr Kruit 
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1988), varies little with position along the major axis. 
This result is independent of galaxy type. 

— The mean levels for the sharpness of the i^T'-band lu- 
minosity peaks indicate that the vertical luminosity 
distributions are more peaked than expected for the 
intermediate sech(z) distribution, but rounder than ex- 
ponential. Since projection of a galaxy causes the pro- 
file to be flattened near the galaxy plane, our result is 
consistent with the hypothesis that all spiral galaxies 
have exponential vertical profiles. The fact that we ob- 
serve this in all our sample galaxies indicates that the 
process at work here is a process intrinsic to the disks 
themselves. 

— For those galaxies with the faintest absolute magni- 
tudes, we find a lack of "rounder" profiles compared 
with those galaxies with brighter absolute magnitudes. 
This may indicate that we are hindered by an under- 
lying dust component, which is concentrated towards 
the galaxy planes, and more extended in the larger 
galaxies than in the smaller ones, at least to an out- 
side observer. It may be that in the smaller galaxies 
this dust component affects relatively fewer data points 
than in the larger ones, thus causing a bias towards 
more sharply peaked vertical profiles. 
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